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ABSTRACT 



We present CO (3-2) interferometric observations of the central region of the 
Seyfert 2 galaxy NGC 1068 using the Submillimeter Array, together with CO (1-0) 
data taken with the Owens Valley Radio Observatory Millimeter Array. Both the 
CO(3-2) and CO(l-O) emission hnes are mainly distributed within ~ 5 arcsec 
of the nucleus and along the spiral arms, but the intensity distributions show 
differences; the CO (3-2) map peaks in the nucleus, while the CO (1-0) emission 
is mainly located along the spiral arms. The CO(3-2)/CO(l-0) ratio is about 
3.1 in the nucleus, which is four times as large as the average line ratio in the 
spiral arms, suggesting that the molecular gas there must be affected by the 
radiation arising from the AGN. On the other hand, the line ratios in the spiral 
arms vary over a wide range from 0.24 to 2.34 with a average value around 0.75, 
which is similar to the line ratios of star-formation regions, indicating that the 
molecular gas is affected by star formation. Besides, we see a tight correlation 
between CO(3-2)/(l-0) ratios in the spiral arms and star formation rate surface 
densities derived from Spitzer 8 fim dust flux densities. We also compare the 
CO(3-2)/(l-0) ratio and the star formation rate at different positions within the 
spiral arms; both are found to decrease as the radius from the nucleus increases. 

Subject headings: galaxies: active, galaxies: individual (NGC 1068), galaxies: ISM, 
galaxies: nuclei, galaxies: Seyfert 
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INTRODUCTION 



The properties of circumnuclear molecular gas (CMG) near an active galactic nucleus 
(AGN) can be influenced by the AGN itself. Some CMG observations of Seyfert galaxies 
with a spatial resolution around 100 to 300 pc have shown that the intensities of C0(2-l) 
line emission a re about two times higher in temperature units than those of CO (1-0) 



emission (e.g., 



Matsushita et al. 



2004 



Hsieh et al. 



20081 ). This is different from the 



properties of molecular clouds in normal or star-forming galaxies, which usually show 
brigh tness temperatures in higher- J transitions simila r to or lower than those in lower- J lin es 



Devereux et al. 



Matsushita et al. 



994 



Mauersberger et al. 



1999 



Dumke et al. 



2001 



Oka et al. 



2007 



20091 ). From the LVG model, the high values of CO(2-1)/CO(1-0) arise 
when the emitting molecular gas is in a low opacity environment and the resulting molecular 
densities are larger than the critical densities of CO(l-O) and C0(2-l) transitions. A 
possible origin for such unusual behavic 
regions due to the strong radiation (e.g., 



2004 



Meijerink &: Spaans 



2005 



r is that CMG near 


AGNs traces X-ra's 


T dominated 


Malonev et al. 


1996; 


Kohno et al. 


2001 


Usero et al. 



Meijerink et al. 



20071 ). Unlike far-UV photon s, X-ray 



photo ns have greater penetration lengths and are more efficient in gas heating ( lUsero et al. 



20041 ). On the other hand, CMG can also b e affected by mechanical processes, such as gas 



entrainment by jets (IMatsushita et al. 



20071). 



Star formation in galaxies is w ell correlated with gas sur: 



the famous Schmidt-Kennicutt law (Schmidt 



1959 



Kennicutt 



ace d ensity as embodied in 



19981 ). In the circumnuclear 



regions (CNRs) and inner structures of galaxies, the molecular gas usually has relatively 
higher densities and temperatures than the molecular gas in galaxy disks; in such 



environme n ts, star formation is expected to be much more vigorous (iKennicutt 



Tan 



2000|). 



1998 



Krumholz fc ThompsonI (120071 ) suggest that the Schmidt-Kennicutt law might 



change slope when the averaged gas density is close to the line critical density. We note that 



-4- 



the critical densities of CO(l-O) (~ 10^ cm ^) and CO(3-2) ~ 3 x 10^ cm ^ are comparable 



with the densities of the CNRs o 
10^ cm"^ 



le.g., 



Matsushita et al 



sorn e nearby active galaxies, which are are around 300 to 



2004 



Sakamoto et al 



2007 



Perez-Beaupuits et al. 



20091 ) 



It is unclear whether or how the Schmidt-Kennicutt law would vary in the CNRs of these 
galaxies when probed with CO(l-O) and CO(3-2) emission. 

Among nearby g alaxies, NGC 1068 is the best-studied prototypical Seyfert. It has a 



distance of 14.4 Mpc flTuUv 



catalog fide Vaucouleurs et al 



19881. : 1" ~ 70 pc) and is classified as (R)SA(rs)b in the RC3 



199ll ). with an inner bar in the central kpc and two tight 



wound spiral arms starting from the tip of the ba r (IScoville et a 



19891 ). The inner 2 kpc is rich in star formation (ITelesco et al. 



1988 



star forming 



Davies et al 



regions are concentrated al ong the spiral arms (e.g. 



1998 



Emsellem 



arms (e.g. 



Myers &: Scoville 



et al 



1987 



Thronson et al 



984), and the active 



Telesco fc Decher 



1988 



2009). Molecu 



Planesas et al 



ar gas is also abundant a 



1991 



Kaneko et al 



1992 



ong the spiral 



Heifer fc Bhtz 



9951) . and a weak o f fset ridge of emission a long the leading side of the bar is also seen 



(iHelfer fc Blitz 



1995 



Schinnerer et al. 



2000l ). which is a typical molecular gas distribution 



in barred galaxies. 



The nucleus of ^ 



Gall 



more et al. 



1996 



GC 1068 shows radio jets (e.g., 



Wilson & Ulvestad 


1983 






Pogge 


1988 




Macchetto et al. 



2004[ 1 and an ionization cone (e.g. 
19941 ). The nuclear optical spectrum has type 2 characteristics, but the polarized spectrum 
shows type 1 features, refle cting the existence of an optically thick torus or disk around the 



central massive black hole (lAntonucci fc Miller 



19851). I ndeed, a sub-parsec scale ion ized gas 



disk perpendicular to the rad io iets has been obser ved (IGallimore et al. 



molecular gas (maser) torus (IGreenhill et al.lll996l ) or disk ( 



19971. 



Gallimore et al. 



20041) ■ with a 



20011) located 



outside of the ioni zed gas disk. Th is maser torus/disk may be surrounded by a pc-scale 



warm dust torus (jJaffe et al 



20041 ). Outside this structure, the existence of a warped 



- 5 - 



mole c ular gas disk is sugges ted by interfere metric CO (2-1) observations (IBaker fc Scoville 



1998 



Schinnerer et al 



20001 ) . but warm molecular gas kinematics dominated by irregular 



(infalling) motions have bee n observed in 2.12 fim 1-0 S(l) molecular hydrogen emission 



(IMiiller Sanchez et al. 



20091 ). Interferometric observations of molecular gas emission with 



angular reso 



(IKrips et al. 



ution s at 0."5-2" also indicated non-circular motions in the central ~100 pc 



20111 ). The intensity ratios of various mole cular species exhibit peculiar values , 



including very high HCj N 



1994 



Kohno et al. 



2001 



/CO and H 



Usero et al. 



CN/HCO+ ratios 



2004 



Krips et al 



Jackson et al. 



2008 



1993 



Tacconi et a. 



Perez-Beaupuits et al 



20091), 



suggesting that the molecular gas 



X-ray emission (lUsero et al. 

2om . 



2004 



very close to t 



Kohno 



2005 



le Seyfert 2 nucleu s is irradiated by stro ng 



Kohno et al. 



2008 



Garcia-Burillo et al. 



The central region of NGC 1068 has been mapped in different CO tranistion lines by 
bo th single-dish a nd int erferometric observati ons. Among t h e form er, it has been observed 



Scoville et al. 



Planesas et al 



'CRAO) 



1989 



Kaneko et al. 



1989 



19831 . CO(l-O), FCRAO)) 

, CO f 1-0) and C0(2- l), IRAM), lYoimg et all 



2009 



Papadopoulos fc Seaquist 



(1999 



C O (1-0 ), Nobeyama) 
(Il99a CO(l-O), 



CO (2-1) and CO (3-2), JCMT), and 



Israel 



CO(l-O), C0(2-l) 



it has been obser ved by 



CO (1-0), NMA), 



, CO (3-2), and CO (4 



Planesas et al. 



Heifer fc Bhtz 



and CO (2-1), IRAM)^ 



Heifer et al. 



(11991 



3), IRAM and JCM' 
CO (1-0), OVRO), 



995LCO(1-0), BIMA), 



(12003 



') . With interferome ters 



Kaneko et al. 



Schinnerer et al 



Krips et al 



(1992 



2000L CO(l-O) 



(12011 



C0(2-l) 



CO(l-O), BIMA), and 

and CO(3-2), SMA/PdBI). In this paper, we present interferometric observations of the 
CO(3-2) line in the central 2 kpc of NGC 1068 using the Submillimeter Array (SMA). 
We also show interferometric observations of the CO (1-0) transition, with uv sampling 
similar to that of the CO (3-2) data, from the Owens Valley Radio Observatory (OVRO) 
Millimeter Array. We describe our observations in Sect. |2]and show the overall molecular 
gas distributions in Sect. [31 Line ratios, molecular gas masses, and kinematics are presented 
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and discussed in Sect. 14. ![ \A.2\ and I4.3[ respectively. The relation between line ratios and 
star formation is discussed in Sect. 14. 4[ and a summary is in Sect. O 



OBSERVATIONS AND DATA REDUCTION 



2.1. SMA CO(3-2) Observations 



We used the SMA to acquire three CO(3-2) datasets on August 13, 15, and 23, 2005. 
The zenith opacity at 225 GHz was about 0.07 in August 13 and 15 and 0.06 in August 
23. Six out of eight 6-m antennas were used in the compact configuration. The receivers 
were tuned to the redshifted CO (3-2) line (344.493 GHz). Correlators were set to cover a 
velocity range of ~ 1700 km s^^ (~ 2 GHz bandwidth) and configured to have a velocity 
resolution of ~ 0.7 km s"^ (= 0.8125 MHz frequency resolution). The SMA antenna 
primary beam has a half power beam width (HPBW) of ~ 36" (~ 2.5 kpc in NGC 1068) 
at 345 GHz. A 7-pointing mosaic (pointings separated by 18" and hexagonal in shape) was 
observed in order to image the inner 72" (5 kpc) of NGC 1068. The phase center was set 
at R.A.= 2M2m40.798s and decl.= -0rf00m47.938s (J2000), corresponding to the AGN 



position (IMuxIow et al. 



19961). 



We performed the data reduction following the standard processes outlined in the 



SMA cookbook^ 



(IScoville et al. 



The visibilities were first calibrated with the IDL-based MIR package 



19931 ) as modified for the SMA. We used 0423 — 013 as the amplitude and 
phase calibrator to track phase and gain variations and used 3C 454.3 for bandpass and 
flux calibrations. 



We calibrated our observations by linearly interpolating the flux densities of the quasar 



^http ://cf a-www . harvard . edu/~ cqi/mircook . html 
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3C 454.3, which were 26.48 ± 1.33 Jy on August 12 and 18.90 ± 0.95 Jy on August 25, 
according to the SMA cahbrator hscl. The expected flux densities of 3C 454.3 were about 
25 ± 1.3 (August 13-15), and 20 ± 1.0 (August 23) during the observations. Mars and 
Uranus were also observed during the observations. However, Mars was resolved by the 
SMA observations, and Uranus was only observed on August 13 and 15. After performing 
the flux calibration with 3C 454.3, we compared the fluxes of Uranus and Mars to check the 
accuracy of our flux calibration. The uncertainty in our flux scale is estimated to be ~ 10%. 

We combined the visibilities from all the mosaic observations and applied the mapping 
task INVERT in MIRIAD to produce "dirty" images with a velocity resolution of 10 km s~^. 
Primary beam correction was taken into account in the mosaic mode of the INVERT 
process. We used natural weighting in the mapping process in order to have the best 
sensitivity. We deconvolved the dirty images with the MOSSDI and MOSMEM packages 
in MIRIAD, and produced moment maps and spectra for further analysis. The resulting 
spatial resolution was 2'.'29 x 2'.'00 (~ 160 x 140 pc^) with a position angle of — 88?8. 



2.2. OVRO CO(l-O) Observations 



We observed the CO(l-O) transition in NGC 1068 (adopted pointing center: aj2ooo 



02:42:40.7 and ,5.i onnn = -00:00:47.7 ) 



millimeter array (jPadin et al. 



1991 



between Apri 


and 


Scott & Finch 


1993) 



and September 1995 using the OVRO 



19931 ). The array comprised six 10.4m 



antennas, which during our observations were deployed in three configuratio ns providing 



a tot al of 40 distinct baselines. We configured the array's digital correlator (iPadin et al. 



19931 ) to provide 112 contiguous frequency channels, each Hanning-smoothed to 4 MHz 
(10.4 km s~^) resolution. We calibrated the data within the OVRO millimeter array 



http : //smal . sma . hawaii . edu/ callist/ callist . html 
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database using the MMA package (IScoville et al.lll993l ). Paired integrations on J0339— 017, 
interleaved with observations of NGC 1068 every 30-40 minutes, were used to correct for 
phase and gain variations; for passband cahbration, we used 3C273, 3C454.3, and/or 
0528+134. Our flux scale was determined by comparing J0339— 017 with Uranus or Neptu ne 



in light of standard bri ghtness temperature models for the latter ( iMuehleman fc Berge 



1991 



Orton et al. 



19861 ). or by bootstrapping from archived observations of the planets (and 



bright, frequently-observed quasars) obtained with similar elevations and coherences. From 
repeatibility of flux measurements on 1-2 week timescales, we estimate that the uncertainty 
in our flux scale is ~ 10%. 



19971 ). we were left with 



After editing for quality in the Difmap package (iShepherdl 
11.4 hours of on-source data. To eliminate contamination by 3 mm continuum emission 
associated with NGC 1068's jet, we subtracted a -ut^-plane model based on the line-free 
channels at the ends of the recorded bandwidth. The line-free uv data were then mapped 
using the IMAGR task in the NRAO AIPS package, with moderately robust weighting 
giving a synthesized beam of 3'.'46 x 2'.'56 (242 pc x 179 pc) at PA 67?5. For deconvolution, 
we adopted a single clean box for all channels that enclosed all of the emission in the 
zeroth moment map. A slightly different reduction of these data was first discussed by 



Baker fc Scovilld (119981 ). Before analyzing moment maps and spectra, we applied a primary 



beam correction appropriate for the ~ 56" OVRO HPBW. 



3. CO Distribution 

In Figure [H we show the integrated intensity and velocity maps of the CO(3-2) line for 
the central region of NGC 1068. CO(3-2) emission appears in the nucleus and along the two 
spiral arms. The CO (3-2) distribution in the nucleus is elongated in an east-west direction, 
and the strongest peak in our data is located ~ 1" east of the nucleus. The distribution 
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of th e CO (3-2) emission i s roughly consistent with that of previous CO (1-0) observations 



[e.g., 



Heifer fc Bhtz 



1995 



Planesas et al 



1991 



Kaneko et al 



1992 



Schinnerer et al. 



20001 ) 



except that the nucleus is relatively brighter in the CO (3-2) image. 

The CO(3-2) intensity distributions in the two spiral arms are different. The southern 
arm is brighter than the northern arm. The second strongest peak in our map is located 
~ 12" south of the nucleus, which is on the southern arm. The locations of these two 
strongest peaks in the map can explain the position of th e peak intensity in an early 



CO (3-2) single-dish image ( jPapadopoulos fc Seaquistlll999l ). which showed a peak intensity 



offset south of the nucleus. The missing flux in the central 14" region is about 20% when 



comp aring our image with a JCMT observation of the central region of NGC 1068 (jlsrael 



20091 ). On the other hand, if we co mpare our data with an area-av eraged CO (3-2) spectrum 



over a circular area of 30" radius (jPapadopoulos fc Seaquist 



19991 ). we find that our map 



only reco vers ~ 50% of the total fl u x. T he missing flux might seem surprisingly large; 



however. 



Papadopoulos fc Seaquist 



fll999l ) have showed that the total CO (3-2) flux is 



actually dominated by a warm diffuse gas phase, which is highly excited and not virialized. 



Besides, similar extended structures of CO (3 



regions of other nearby AGNs, such as M51 ( iMatsushita et al. 



2) emission are also o 



Dserved in the central 



200J). Most of this diffuse 



gas is likely resolved out and not detectable in our interferometer map. 

The CO(l-O) integrated intensity and velocity maps are shown in Figure [2l As in the 
CO(3-2) map, CO(l-O) emission is located in the nucleus and along the two spiral arms, but 
the strongest peak is located at the southern arm, and the peak intensity at the core region 



is relatively weak. Our CO (1-0) image is in 



CO (1-0) observat ions (jPlanesas et al. 



Heifer et al. 



1991 



general agreement w i th previous inter: 



Kaneko et al. 



1992 



Schinnerer et al. 



erometric 



2000 



20031 ). However, the BIMA image seems to show more emission at interarm 



regions. The difference between our results and the BIMA image might be due to the fact 
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RA (J2000) RA (J2000) 

Fig. 1. — Left: CO (3-2) integrated intensity (moment 0) map of the central region of 
NGC 1068. The grayscale range is shown in the wedge at right in units of Jy beam~^ km s^^. 
The contour levels are 10, 20, 30, 40, 50, 70, 90, 110, and 130 xa, where la = 
4.0 Jy beam"^ km s"^ The synthesized beam is 2'.'29 x 2'.'00 (~ 160 x 140 pc^) with a 
P.A. of — 88?0, which is shown in the bottom left corner. The cross indicates the galacti c 



center determined from the 5 GHz and 22 GHz radio continuum data (iMuxlow et al.lll996h . 
Right: CO (3-2) intensity- weighted velocity field (moment 1) map. The grayscale range is 
shown in the wedge at right from -180 to 180 km s"^. The contour levels range from —150 
to +150 km s"^ with a 30 km s~^ interval. 



- 11 - 




RA (J2000) HA (J2000) 



Fig. 2. — Left: CO (1-0) integrated intensity (moment 0) map of the central region of 
NGC 1068. The grayscale range is shown in the wedge at right in units of Jy beam~^ km s^^. 
The contour levels are multiplies of 5a, where la = 0.81 Jy beam~^ km s~^. The synthesized 
beam is 3'.'46 x 2'.'56 with a P.A. of 67?5, which is shown in the bottom left corner. The cross 
is the same as in Fig. [2j Right: CO(l-O) intensity weighted velocity field (moment 1) map of 
the central region of NGC 1068. The grayscale range is shown in the wedge at the right, from 
-180 to 180 km s~^. The contour levels range from —150 to +150 km s^^ with a 30 km s~^ 
interval. 
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that the BIMA image have included modelled visibilities derived from single-dish NRAO 
observations and might contain large-scale emission that would be resolved out by our 
interferometric observations. 



The corresponding m issing flux o: 



with the BIMA results of 



Heifer et al. 



CO(l -O) in the core region is about 20% comparing 



( l2003l ). For the spiral regions, the situation is more 



complicated. We have estimated the missing fluxes for several selected spiral regions (see 
Figure 3) and obtained diverse results. For example, the missing flux for R4 is about 70% 
and for R17 is about 38%; however, we detect more flux than BIMA for R15, so there should 
be no missing flux for R15. We note that R4 is an inter-arm region, which is dominated by 
diffuse emission, while R15 is around a compact structure peak and R17 is at the boundary 
of an arm structure. These different results might be caused by different uv coverage of the 
BIMA and our observations. 

The overall CO (1-0) velocity fleld is similar to that of CO (3-2). Both have a rotation 
axis with P.A. about —30° and a rotation velocity of 120 km s~^ in the nucleus. In the 
spiral arms, both kinematic major axes run from east to west with a rotation velocity of 
180 km s-^ 



4. Discussion 

4.1. CO Line Ratios 

Before deriving line ratios at various regions, we flrst matched the uv range between our 
CO(3-2) and CO(l-O) datasets. In Figure [3l we overlay the CO(3-2) (solid contours) and 
CO(l-O) (dashed contours) intensity distributions. In this image, the shortest uv length is 
set to 16 kA, and the image resolutions are convolved into the same resolution (3'.'46 x 2'.'56 
with a P.A. of 67?5). As mentioned above, the overall distributions of both CO (3-2) and 
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CO(l-O) are very similar, but some of the intensity peaks have shifted positions with 
respect to each other. The most obvious example is in the southern spiral arm, about 
12" — 14" south of the nucleus, where the CO(3-2) peak is located in the inner part of the 
spiral arm, but the CO(l-O) peak is shifted toward the outer edge of the spiral arm. 

We divide the central region of NGC 1068 into 25 regions covering the nucleus 
and spiral arms, as shown in Figure [3] (CI and Rl - R24). The size of each region is 
4" X 5" except CI, which is 6" x 5". The intensity scale is then converted to a brightness 
temperature scale, and the CO (3-2) /CO (1-0) brightness temperature ratios, R31, are 
derived for each region using the MIRIAD task MATHS. Area- averaged spectral peak and 
integrated brightness temperatures and R31 for each region are shown in Table [H 

The central core region CI exhibits a very high integrated intensity ratio with 
i?3i = 3.12 ± 0.03 and a spectra peak ratio of 2.83 ± 0.10. Figure HI sho ws the spectra of C I. 



The CO(3-2)/CO(l-0) line ratio is slightly smaller than the re s ults o f 



Krips et al. 



Krips et al 



( 1201 ih 



(120111 ) have a much higher 



which are about 4-6. However, the observations of 
angular resolution and might contain relatively more contribution from the circumnuclear 
region of the AGN, while our result is an average value over a much larger area and are 
likely to include emission from outside regions. In fact, if we only consider the line ratio 
within the brightest be am at the ce r iter, w e would obtain a line ratio of ~ 4.6, which is 



similar to the results of 



Krips et al 



(1201 ih 



On the other hand, the spiral arm regions Rl - R24 have a wide range of integrated 
intensity ratios of 0.24 - 2.34 and an average value of 0.75 with standard deviation 0.47. 
We note that the wide range of the intensity ratios might be caused by the different spatial 
distributions of the CO(l-O) and CO(3-2) emission due to varying physical conditions of 
the molecular gas. There is no obvious ratio difference between the northern and southern 
arms; the average ratio of the northern arm regions (Rl - R4, R18 - R24) is 0.8 and that 
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o t rr 



-00 0106 - 



2^42™42f0 4lf5 



40fO 



39?5 



39?0 



4lfO 40f5 
RA (J2000) 

Fig. 3. — Integrated intensity maps of CO (3-2) (solid contours) and CO (1-0) (dashed 
contours), overlaid on the continuum-subtracted HST F658N image (greyscale) of the central 
region of NGC 1068. Solid contour levels for CO(3-2) are 10, 30, 50, 70, 90, 110, and 
130 X 5.3 Jy beam"^ km s~^, and dashed contour levels for CO(l-O) are 2, 4, 6, 8 x 
5 Jy beam"^ km s~^. The CO(3-2) and CO(l-O) data are matched to the same uv range 
and have the same synthesized beam size of 3'.'46 x 2'.'56 (~ 242 x 179 pc^) with a P.A. of 
67?5, which is shown in the bottom-left corner of the image. Cross is the same as in Figure [H 
We also plot 25 boxes (CI and Rl - R24) that are used to calculate the line ratios. 
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of the southern arm regions (R5 - R17) is 0.7. However, there is a difference between inner 
arm regions and outer arm regions; the average ratio of the inner arm regions (Rl - R4, R9 
- R17, and R19 - R22) is 0.88 ± 0.14, while the outer arm regions (R5 - R8, R18, R23, 
and R24) show 0.42 ± 0.05. This difference indicates that there is a large-scale gradual 
decrease of the line intensity ratio from the inner radii to the outer radii. The radial 
change of molecular line ratios, and therefore of the physical con ditions of the molecular 



gas^ 



has also been detected in the Milky Way and other galaxie s ( 



1994 



Sakamoto 



Matsushita et al 



994 



Sakamoto et al 



1997 



Petitpas &: Wilson 



Turner 



20001 : 



1993 



Aalto et al. 



Paglione et al. 



2001 



20101), suggesting that this trend is common. 



The line ratios -R31 in the spiral arm r egions are similar to th ose observed in the centers 



of nearby normal and starburst galaxies. 



Devereux et al. 



(119941 ) observed seven nearby 



starburst galaxies and found that the ratios are in the range of 0.5 - 1.4 with an average 



ratio of 0.64 ± 0.06. A survey toward the centers of 28 ne arby star-forming ga . 



axies showed 



19991 ). On 



that most of their R31 are within the range of 0.2 - 0.7 (iMauersberger et al. 
the other hand, the value of R31 in the nucleus is much larger than those observed in the 
nuclei of the nearby galaxies mentioned above, albeit similar to that observed in the nucleus 
of the Seyfert 2 galaxy M51. The inner -uf -truncated -R31 of the central molecular core of 
M51 also has a very high value of ~ 5.2 ± 1.7 wi thin a beam size of 3^.^9 x 2'.'6 or 160 pc x 
110 pc (missing flux corrected R31 is > 1.9 ±0.2; 



Matsushita et al 



20041 ). If we only correct 



the 20% missing flux of the CO(l-O) emission in the cor e region of NGC 1068 an d ignore 



Matsushita et al. 



(120041 ) ■ the R31 



the missing flux of the CO (3-2) following the method of 
ratio of the core region will become 2.50 instead of 3.12. This suggests that the physical 
conditions of molecular gas around the Seyfert 2 AGN in NGC 1068 are very different from 
those in the centers of nearby star-forming and starburst galaxies, but might be similar to 
those around the Seyfert 2 nucleus of M51. 
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Fig. 4. — Area-averaged spectra within the CI region. The sohd hne shows the spectrum of 
CO(3-2), and the dash hne shows the spectrum of CO(l-O). The uncertainties per channel 
are about 0.01 K for the CO(3-2) spectrum and 0.02 K for the CO(l-O) one. 



Table 1. Peak and integrated brightness temperature of CO (3-2) and CO (1-0), and the 
CO(3-2)/CO(l-0) hne ratios, in different areas. 



Area Peak CO (3-2) Peak CO (1-0) Peak R31 f TsiCO J = 3 - 2) dv f TsiCO J = 1 - 0) dv R31 X offset Y offset 

[K] [K] [K km s~^] [K km s~^] [arcsec] [arcsec] 



CI 


1.73 ±0.01 


0.61 ± 0.03 


2.83 ±0.10 


438.5 ± 0.7 


140.6 ± 1.5 


3.12 ± 0.03 


-1 





Rl 


1.80 ± 0.03 


1.55 ± 0.06 


1.16 ± 0.05 


103.76 ± 0.97 


44.3 ± 1.8 


2.34 ± 0.10 





14 


R2 


1.53 ± 0.05 


1.90 ± 0.08 


0.80 ± 0.04 


104.97 ± 1.54 


133.4 ± 2.8 


0.79 ± 0.02 


4 


14 


R3 


1.29 ± 0.06 


1.98 ± 0.05 


0.65 ± 0.03 


100.06 ± 2.07 


149.8 ± 1.8 


0.67 ± 0.02 


8 


12 


R4 


0.43 ± 0.03 


0.73 ± 0.04 


0.58 ± 0.05 


6.26 ± 0.81 


26.6 ± 1.2 


0.24 ± 0.03 


12 


9 


R5 


0.78 ± 0.04 


1.92 ± 0.11 


0.41 ± 0.03 


56.74 ± 1.39 


90.0 ± 3.8 


0.63 ± 0.03 


16 


6 


R6 


1.54 ± 0.05 


3.65 ± 0.15 


0.42 ± 0.02 


78.95 ± 1.45 


131.1 ± 4.4 


0.60 ± 0.02 


18 


9 


R7 


0.71 ± 0.04 


1.97 ±0.14 


0.36 ± 0.03 


25.72 ± 1.23 


93.0 ± 4.4 


0.28 ± 0.02 


18 


14 


R8 


0.84 ± 0.03 


1.80 ± 0.08 


0.47 ± 0.03 


28.85 ± 1.09 


94.4 ± 2.5 


0.31 ± 0.01 


20 


3 


R9 


1.35 ± 0.03 


3.63 ± 0.08 


0.37 ± 0.01 


108.83 ± 1.11 


199.9 ± 2.9 


0.54 ± 0.01 


16 


-1 


RIO 


1.49 ±0.07 


2.20 ± 0.08 


0.68 ± 0.04 


62.03 ± 1.92 


111.0 ±2.2 


0.56 ± 0.02 


12 


-5 


Rll 


1.49 ±0.07 


1.94 ± 0.03 


0.77 ± 0.04 


94.26 ± 2.63 


172.2 ± 1.1 


0.55 ± 0.02 


8 


-8 


R12 


1.74 ± 0.06 


3.02 ± 0.04 


0.57 ±0.02 


112.76 ±2.01 


176.3 ± 1.4 


0.64 ±0.01 


4 


-10 


R13 


1.93 ± 0.09 


1.66 ± 0.07 


1.16 ±0.07 


194.32 ± 3.64 


139.4 ± 2.8 


1.39 ± 0.04 





-12 


R14 


2.46 ± 0.09 


3.40 ± 0.03 


0.72 ± 0.03 


252.08 ±3.58 


192.5 ± 1.3 


1.31 ±0.02 


-4 


-13 


R15 


2.68 ±0.08 


3.18 ± 0.09 


0.84 ±0.04 


239.52 ±3.24 


374.4 ±3.4 


0.64 ±0.01 


-8 


-11 


R16 


3.46 ±0.11 


4.80 ±0.18 


0.72 ± 0.04 


211.39 ±3.69 


301.1 ±6.0 


0.70 ± 0.02 


-12 


-8 


R17 


1.67 ±0.01 


1.68 ± 0.06 


0.99 ± 0.04 


106.98 ± 0.22 


141.2 ±2.3 


0.76 ± 0.01 


-16 


-5 


R18 


0.79 ±0.01 


2.47 ±0.07 


0.32 ±0.01 


56.68 ± 0.62 


159.4 ±2.8 


0.36 ± 0.01 


-20 


1 


R19 


1.28 ±0.05 


1.72 ± 0.08 


0.75 ± 0.04 


49.41 ± 1.41 


75.5 ± 2.4 


0.65 ± 0.03 


-16 


6 


R20 


1.94 ± 0.07 


3.19 ± 0.06 


0.61 ±0.02 


115.44 ± 2.06 


168.9 ± 1.8 


0.68 ±0.01 


-12 


10 


R21 


1.83 ±0.05 


1.66 ± 0.06 


1.10 ±0.05 


97.04 ± 1.48 


81.3 ± 1.8 


1.19 ±0.03 


-8 


13 


R22 


1.10 ±0.04 


0.47 ± 0.04 


2.32 ±0.22 


62.59 ± 1.26 


46.3 ± 1.3 


1.35 ± 0.05 


-4 


14 
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4.2. Molecular Gas Column Density and Mass 



Because CO(l-O) traces the bulk of the molecular gas, we can calculate the molecular 
gas column density and mass from the CO (1-0) integrated intensities. The column density 
of the molecular hydrogen can be estimated with a conversion factor, Xqq, 



iVn, = Xco y Tb(CO J = 1 - 0) dv [cm~ 



wher e Xcn = 0-4 x 10^° cm ^ ( K km s ^) ^ for circumnuclear molec ular gas (jWilson 



1995 



Mauersberger et al. 



1996 



WeiB et al. 



2001 



Espada et al. 



2009h and Xc.n 



1987 



Young et al. 



199l|). 



3 X 10^° cm~^ (K km s~^)~^ for spiral arm regions (ISolomon et al. 
The column density A^H2 for each region is calculated using J Tb(CO J = 1 — 0) dv shown 
in Table [H and the derived values are shown in Table |2l We also derive the molecular gas 



mass, , which is displayed in Table |2i The total missing flux of C0( 



central circular area of 30" radius is about 35% (jPapadopoulos fc Seaquist 



-0) w ithin the 



19991). However, 



the missing flux is dominated by extended structures, so the true missing fluxes in the 
core and the compact spiral structures should be smaller. For example, the missing flux of 
CO (1-0) in the core region is about 20% as discussed in Sect. O We thus expect that the 
errors of the column densities caused by the missing flux are less than 35%. 



Table 1 — Continued 



Area 


Peak CO{3-2) 


Peak CO(l-O) 


Peak R31 


/Tb(CO J = 3-2) dv 


/Tb(CO J=1-0) dv 


R31 


X offset 


Y offset 






[k: 




[K km s 1] 


K km s ''-] 




[arcsec] 


[arcsec] 


R23 


1.02 ±0.05 


2.72 ± 0.07 


0.38 ± 0.02 


45.12 ± 1.70 


94.9 ± 2.1 


0.48 ± 0.02 


-24 


-4 


R24 


0.96 ±0.03 


2.95 ±0.16 


0.33 ±0.02 


39.06 ± 0.86 


126.5 ±4.5 


0.31 ± 0.01 


-24 


-9 



Note. — Column (1): Region name. Column (2): Peak brightness temperature of CO(3— 2). Column (3): Peak brightness 
temperature of CO(l-O). Column (4): Line ratio of peak brightness temperature of CO(3-2)/CO(l-0). Column (5): Integrated 
brightness temperature of CO(3-2). Column (6): Integrated brightness temperature of CO(l-O). Column (7): Line ratio of 
integrated brightness temperature of CO(3-2)/CO(l-0). The errors Column (8,9): Position offset of region center from the 

phase center. Errors were estimated from the statistical errors of the spectra; uncertainties of our flux calibration (about 10% | 

I L 

for both CO(l-O) and CO(3-2) lines) and missing fluxes are not included. CO 

I 
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Table 2. Molecular gas column density and mass, 8 fxm. dust and Ha + [N II] emission, 

and star formation rate surface density. 



Area 






/s^m.dust 


/Ha + [NII] 


SsFR 




Xl02l [cm-2] 


Xl06 [Mq] 


[MJy sr-l] 


xlO"* [eps pixel"^]'' 


[Mq yr-1 kpc-2] 


CI 


5.6 ±0.3 


12.3 ±0.8 


1.05 X 10* ± 6.4 


73.5 


62.4 ± 0.04 


Rl 


13.3 ±0.5 


19.5 ± 0.8 


2,38 ± 6.4 


1.51 


1.41 ±0.04 


R2 


40.0 ± 0.8 


58.6 ± 1.2 


260 ± 6.4 


> 1.87'' 


1.54 ±0.04 


R3 


44.9 ± 0.5 


65.8 ±0.8 


337 ± 6.4 


> 2.57'' 


2.00 ± 0.04 


R4 


8.0 ±0.4 


11.7 ±0.5 


241 ± 6.4 


2.35 


1.43 ± 0.04 


R5 


27.0 ± 1.1 


39.5 ±1.7 


180 ± 6.4 


1.08 


1.07 ± 0.04 


R6 


39.3 ± 1.3 


57.6 ±1.9 


166 ± 6.4 


> 0.81'' 


0.98 ± 0.04 


R7 


27.9 ± 1.3 


40.9 ± 1.9 


107 ±6.4 


b 


0.63 ± 0.04 


R8 


28.3 ±0.8 


41.5 ± 1.1 


89 ± 6.4 


> 0.72'' 


0.52 ± 0.04 


R9 


60.0 ±0.9 


87.9 ± 1.3 


163 ± 6.4 


> 0.65'' 


0.97 ±0.04 


RIO 


33.3 ±0.7 


48.8 ± 1.0 


158 ± 6.4 


0.77 


0.94 ± 0.04 


Rll 


51.7 ±0.3 


75.7 ±0.5 


174 ± 6.4 


0.89 


1.03 ± 0.04 


R12 


52.9 ±0.4 


77.5 ±0.6 


226 ± 6.4 


0.97 


1.34 ±0.04 


R13 


A'\ o _l_ n o 
41.5 ± 0.8 


61.3 ± 1.2 


212 ± 6.4 


0.76 


l.ZD ± 0.U4 


R14 


57.8 ± 0.4 


84.6 ±0.6 


236 ± 6.4 


1.21 


1.40 ±0.04 


R15 


112.3 ±1.0 


164.5 ± 1.5 


311 ±6.4 


2.17 


1.84 ± 0.04 


R16 


90.3 ± 1.8 


132.3 ± 2.6 


381 ± 6.4 


3.49 


2.26 ± 0.04 


R17 


42.4 ±0.7 


62.1 ± 1.0 


345 ± 6.4 


2.81 


2.05 ± 0.04 


R18 


47.8 ±0.8 


70.0 ± 1.2 


172 ± 6.4 


0.62 


1.02 ± 0.04 


R19 


22.7 ±0.7 


33.2 ± 1.0 


185 ± 6.4 


1.15 


1.10 ±0.04 


R20 


50.7 ±0.6 


74.2 ± 0.8 


207 ± 6.4 


1.45 


1.23 ±0.04 


R21 


24.4 ±0.5 


35.7 ± 0.8 


260 ± 6.4 


1.33 


1.54 ± 0.04 


R22 


13.9 ±0.4 


20.4 ± 0.6 


231 ±6.4 


1.37 


1.37 ±0.04 


R23 


28.5 ±0.6 


41.7 ±0.9 


123 ± 6.4 


0.27 


0.73 ± 0.04 


R24 


38.0 ± 1.4 


55.6 ± 2.0 


105 ± 6.4 


0.15 


0.62 ± 0.04 



Note. — Column (1): Region names. Column (2): Column density. Column (3): Moleculcir gas mass. Column (4): Average 
intensity of dust emission. Column (5): Average intensity of Ha + [Nil] emission. Column (6): Star formation rate surface 
density. 
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^Image units in electrons per second per pixels (eps pixel ^) 

''These areas do not have complete Ha emission information, since these regions are either located at the edge of or outside 
the HST image. 
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In the spiral arm regions, the box sizes correspond to 280 x 350 pc^. The average 
values of A^H2 ^iid M^^ per region are 39.8 x 10^^ cm^^ and 58.9 x IO^Mq. The standard 
deviations of A^H2 ^ind are large inside the spiral arms, with a value of 24.0 x 10^^ cm~^ 
and 34.8 x lO^M©, respectively, indicating that the molecular gas content within the spiral 
arms varies substantially from region to region. On the other hand, A^H2 and Mn^ for the 
central core CI are estimated to be 5.6 x 10^^ cm~^ and 12.3 x 10^ Mq (note that the area 
of CI is ~1.5 times larger than those in the spiral arms). This indicates that the column 
density and the mass of the central region are much smaller than the average values in the 
spiral arms. However, we note that the R^i of CI is very different from those of the spiral 
arm regions (see Sect. 14. ip, and we have us ed different conversion factors in estimating the 



mass and column density (IWeiB et al. 



200l|). 



Tacconi et al.l (Il994j ) and [Sternberg. Genzel fc Tacconil (Il994j ) showed that the A^H2 



column density in the core region is about 4 x 10^^ cm~^ using the Galactic CO to A^H2 
conversion factor. Our conversion factor for the core region is about | of the Galactic 
conversion factor and we obtained a A'h2 column density of 5.6 x 10^^ cm~^; in other 
words, the different results are mainly caused by the different conversion factors ad opted. 



1995 



Mauersberger et al. 


1996; 


WeiB et al. 


2001; 


Espada et al. 


2009) 



20091 ). In particular, the 



conversion factor in the center of NGC 1068 could be six times lower than the Galactic 



value ( lUsero et al. 



2004 



Garcia-Burillo et al. 



2010h . 



In Figure El we compare our CO(3-2 



Chandra X-ray Observatory (lYoung et al. 



image with the X-ray image obtained with the 



200ll ). The X-ray image displays clear emission 



from the ionization cone emanating from the nucleus of NGC 1068; however, we note that 
there is an obvious dimmed area at the center of the image. This dimmed area matches 
well with the central core of our CO (3-2) image, suggesting that the obscuring material of 
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the X-rays is closely related to the molecular gas. 



4.3. CO (3-2) Rotation Curve 



The rotation curve in the nucleus of NGC 1068 is plotted in Figure El The data points 
and the error bars were derived using the MIRIAD task VELFIT. We ignore points within 
30° of the minor axis ( P. A. = 60°), and we assume the inclination angle to be 45° following 



Schinnerer et al. 



(|2000[ ). Within the central 2" (~ 140 pc), the rotational velocity increases 
with radius (rigid rotation). The enclosed mass within the central 2" can be estimated 
using M(r < 2") ~ rv^ j G, and we find th a t the total mass within this area is 3.7 x 10* M©, 



which is consistent with 



Schinnerer et al 



(I2OOOI ). The molecular gas mass in CI is about 



12.3 X 10^ Mq, so that the gas-to-dynamical mass r atio is 3%. This value is a factor q 



lower than those in star forming galaxies (~ 10%) (iSakamoto et al. 



1999 



Koda et al 



a few 



20051). 



For radii larger than 2", the rotation curve be comes s 



curve is well fitted with the Brandt rotation curve (iBrandt 



owly decreasing. This rotation 



1960): 



V 



R 



1 _|_ 2 / R 



_3_ 
2n 



(2) 



3 3 \ -Rmax 

with Vmax = 116 km s^^, -Rmax = 2'.'7, and n = 6 (Figure [6]). The mass from the Brandt 
model is: ^ 

Mto, = (0 " Yi,^, (3) 

The total mass inside 190 pc (~ 2.7") radius is 7.2 x 10*^ M©. 



4.4. Relation between CO Line Ratio and Star Formation 



As shown in Sect. 14.11 i?3i va ries along the 



closely related to star formation ( iKomugi et al. 



spira. 



20071 ) 



arms. Since the CO (3-2) line is more 
and i?3i increases with increasing 
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0100 
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41?0 40?5 
RA (J2000) 



Fig. 5. — Comparison betwe en our CO (3 



Chandra X-ray Observatory ( jYoung et al 



2) im age and the X-ray image taken with the 



200ll ). The contours show CO (3-2) emission, 



with levels 10, 20, 30, 40, 50, 70, 90, 110, and 130 x 4.41 Jy beam'^ km s-\ The gray scale 
shows the X-ray emission from to 1000 counts pixel~^(0.4 — 5.0 keV). 
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Radius (arcsec) 



Fig. 6. — Rotation curve of CO(3-2) in the central region of NGC 1068. Crosses with error 
bars are the data points, and the sohd curve is the fitted Brandt rotation curve (see the 
main text for details). 
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star formation efficiency (IMuraoka et al.l 120071 ). the variation in R31 along the spiral arms 
of NGC 1068 may also be related to star formation. To study the cause of the variation, we 
compare R^i with star formation surface density along the spiral arms. 

We use both infrared and optical observations as star formation tracers. The infrared 
data are o btaine d from the Spitzer Space Telescope IRAC 3.6 /im and 8 /im imageqf . 



WuetaL 



(I2OO5I ) showed that dust emission at 8 /im can be used as a star formation 
indicator. We estimate the dust emission of in the spiral arms of NGC 1068 using the 
observed Spitzer IRAC 3.6 /im and 8 fim fluxes: 

/8/im(dust) = /s/xm — Vsfa.eum, 



(4) 



where rjs = 0.232 ( iLeitherer et al.lll999l : iHelou et al.l 12004 ). The derived values are shown 
in Table H 

Ha emission is also considered as a tracer of star formation. For NGC 1068, the 
wavelength of Ha is shifted to 6587.8 A and those of [N II] to 6573, 6603A. To obtain 
the Ha emission, we retrieved the Hubble Space Telescope (HST) F658N and F791W 
images from the STScI archive!. The F658N narrow-band ffiter has a central wavelength of 
6590.8 A with a bandwidth of -30 A; therefore, the F658N image includes the Ha line as 
well as [N II] lines, which are shifted to 6573 and 6603Afor NGC 1068. Since this image 
also includes continuum emission, we correct it using the nearby broad-band ffiter F791W. 
The F791W ffiter has a central wavelength of 7881 A with a bandwidth of ~1231 A. The 
continuum subtraction can be performed using the following equation: 

/nb^Abb — /bb^Anb 



/i 



line 



AA 



BB 



AA 



(5) 



NB 



http : / / ssc . spitzer . caltech . edu/ archanaly/ archive . html 



http : / /archive . stsci . edu 
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where /une is the continuum-subtracted line flux, /nb and /bb are the observed total 
fluxes in the narrow-band and broad-band filters, respectively, and AAnb and AAbb are 
the bandwidths of the narrow-band and broad-band filters, respectively. The resulting 
continuum-subtracted F658N image can be considered as the Ha + [N II] image. The 
derived values are shown in Table [2l We cannot separate the Ha from the [N II] lines 
and we cannot correct the Ha + [N II] image for internal extinction with available data. 
Therefore, we can not obtain any quantitive information from this image; we only use the 
continuum-subtracted emission as a possible indicator of relative star formation rates. 

We compare the HST Ha + [N II] line intensity with the Spitzer 8 fim dust intensity. 
As can be seen in Figure [71 the Ha -|- [N II] line intensity and the 8 fim dust intensity 
are linearly correlated, indicating that the former image is effectively also a good star 
formation tracer. We note that the CO distribution also matches very well with the dust 
lanes along the spiral arms in the HST Ha + [N II] line image (Figure [3]). The star forming 
regions are also located along the spiral arms but slightly shifted toward the outside of the 



spiral arms, especially in the southern spiral arm. This is si milar to the resu 



observations for t 



1999 



le sp iral arms of the nearby galaxy M51 (jVogel et al. 



Koda et al.l 



1988 



t s of previou s 



Aalto et al. 



20091). 



We also compare molecular gas properties with the star formation rate surface density 
I^Ssfr) derived from the Spitzer IRAC 8 /im and 3.6 fim data. The Ssfr is derived from 
the 8 fim dust luminosity: 

SsFR _ uL^[8fim{dust)] 



f lWu et al 



{Me yr-i; 



1.57 X 109 L 







20051 ). where Ly is the 8 /xm dust luminosity. The derived values are shown in 



Table H 

Figure M shows the comparison between R^i and Ssfr of the spiral arm regions. The 
correlation coefficient r between the R31 and Ssfr is 0.33 with a probability p = 89% of 
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Fig. 7. — Comparison of the Spitzer 8 /iin dust intensity with the HST continuum-subtracted 
Ho; + [N II] Una intensity. The HST line intensity is not flux-cahbrated and is expressed 
in instrument units. The data points with upper arrows do not have complete Ha emission 
information because they are located either at the edge of or outside the HST image. 
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3 



CO 







0.5 1.0 1.5 2.0 2.5 
SFR(M0yr"^kpc"^) 



Fig. 8. — Comparison between the CO(3-2)/CO(l-0) integrated intensity ratio, R^i, and 
the star formation rate surface density (Ssfr)- The error bars have included the ~ 10% 
uncertainty of our flux calibrations for the CO (3-2) and CO (1-0) emission. 
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Fig. 9. — Relation between the surface density of molecular gas Sgas and Ssfr for normal 
galaxies (open squares), starburst galaxies (open circles), and the spiral arm regions of 
NGC 1068 (filled circles). The gas surface density is derived from CO( l-O) emission. The 
normal and starburst galaxy samples are obtained from iKennicuttI (1l998l ). The solid triangle 
is the nuclear region of NGC 1068 (CI). The solid line represents the power-law fit to all 
data points. 
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Fig. 10. — Relation between integrated CO (3-2) and FIR luminosities. The filled circles 
represent the data of the spiral arm regions of NGC 1068, and the solid tria ngle is the 



MaoetaL 



nuclea r region of NGC 1068 (CI). The open circles are the data obtained from 
(120101 ). The solid line represents the linear fit to all data points. The outlier point with the 
lowest CO (3-2) emission is R4. 
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mutual correlation. We note that most of the Rsi of nearby star- 



0.2-0.7 (IMauersberger et al. 



Muraoka et al. 



brming galaxies are within 



(120071 ) found that the CO 



19991). For examples, 

emission of the starburst galaxy M83 has R31 < 1 and show a good correlation between the 
i?3i and SFE of their data. If we consider the data with R^i < 1 in our results, we find 
that the correlation between the R31 and Ssfr becomes highly significant with r = 0.63 
and p = 99%. This suggests that in normal star-forming regions the physical conditions 
of molecular gas indicated by -R31 are well correlated with dust emission. The data points 
where R31 > 1.0 are obviously out of the correlation, suggesting that the large R31 might 
be caused by other reasons, such as a different heating mec hanism for the molecular gas or 
different distributions of warm and cool molecular gas (e.g., IHo et al.l 119871 ). 



Figure [9] shows the relation between the surface density of molecular gas Sgas and 
SsFR- The star formation rate surface densities of the spiral arm regions of NGC 1068 
are much higher than those of normal galaxies and similar to those of starbur st galaxies; 



19981). This result 



however, they all seem to follow the same star formation law (iKennicuttI 
strongly supports the idea that the spiral arms in the inner ~ 2 kpc region are experiencing 
a starburst. On the other hand, the molecular gas at CI is obviously offset from the 
Schmidt-Kennicutt law, suggesting that CI is mainly affected by AGN activities instead of 
star formation. 

In Figure \TU\ we present the Schmidt-Kennicutt law for the CO (3-2) emission. We 
used the FIR lumi nosity instead of the SFR in Figure |9] so that it is easier to compare 
with the results of 



Mao et al. 



(120 lOf). We first derive the star formation rate from the 



observed 8/im dust luminosity (IWu et al. 
SFR(M0?/r-i) - 1.7 x 10-^°{Lfir/ Lq) JKennicutt 



2005r) and then use the SFR-Lpj^ relation. 



19981), to derive LpjR. We note that 



this is effectively a Schmidt-Kennicutt law for CO (3-2) emission since Lpm is proportional 
to SFR, and the integrated CO (3-2) luminosity represents the molecular mass in relatively 
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warm and /o r dense reg i ons. We find that most of our results follow the same relation 



obtained by 



Mao et al. 



(120 lOf ) for nearby galaxies. The only outlier point, which has the 



lowest CO ( 3-2) e mission, is the interarm region R4. When we combine our results with 



Mao et al. 



( I2OIOI ). we find that the derived power-law index of the Lco(3-2)^-^F/i? relation 
is ~ 0.9. This value is smaller than the index of the traditional Schmidt-Kennicutt law, 
which is a round 1.0 to 2.0, but is reasonable for the excitation conditions of warm and/or 



dense gas (IKrumholz fc Thompson 



2007 



Narayanan et al. 



20081). 



the NGC 1068 data is about 



Kennicutt 



( 11998h and 



Mao et al 



In Figure and Figure [TDl the physical scale of 
280 X 350 pc^, which is much smaller than those of 

(120101 ). We note that there is no correlation between Egas and Ssfr in Figure |9] when only 
considering the NGC 1068 data. On the other hand, the CO (3-2) emission from the same 
NGC 1068 regions show a very good correlation with the infrared as shown in Figure [THl 
the correlation coefficient r is 0.789 with a probability p = 99.9%. The power-law index of 



the L, 



co(3-2)^LpfR relation for the NGC 1068 data alone (excluding R4) is about 0.49. This 
value is significantly lower than a typical power-law index of the Schmidt-Kennicutt law. 

In Figure [T0| the CO (3-2) to FIR luminosity distribution is generally consistent with 
previous studies, which have a power law index of the Schmidt-Kennicutt law ~1.0. On the 
other hand, the Schmidt-Kennicutt law derived from the CO(l-O) data has a steeper slope 
than that from the CO(3-2) data (Figure [9]). Furtheremore, if we only consider our NGC 
10 68 data, we find a very fi at po wer law index, which cannot be explained by the model 



of 



Krumholz &: ThompsonI ( 120071 ) with a different critical density. One possibility is that 



the gas is under sub - therm al conditions, which would produce a flatter KS law as shown 



by 



Narayanan et al. 



( I2OO8I ). This interpretation is also consistent with the relatively large 
line ratios of CO (3-2) /CO (1-0). We also note that most of our data are under the average 
value of the KS law, indicating that the SFE of the inner spiral regions of NGC 1068 is 
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smaller than that in most of the sources in Figure [TOl 

There is also a radial variation of star formation activity in the spiral arm regions. The 
average star formation rate surface density of the inner arm regions (Rl - R4, R9 - R17, 
and R19 - R22) is 1.45 M© yr^^ kpc^^ and the outer arm regions (R5 - R8, R18, R23, and 
R24) is 0.80 Mq yr^^ kpc^^. In other words, the radial variation of the physical conditions 
of the molecular gas mentioned in Sect. 14.11 is correlated with the radial variation of the 
galaxy's star formation. 

5. SUMMARY 

We have shown and compared the emission of different CO rotational transitions of 
the prototypical Seyfert 2 galaxy NGC 1068 observed with millimeter and submillimeter 
interferometers. The molecular gas in the central part of this galaxy is distributed in a 
central core and outer spiral arms. Both the CO (1-0) and CO (3-2) lines show similar 
distribution along the spiral arms, and most of the molecular gas mass is located in the 
spiral arms. However, the nucleus is rather different; the strongest CO (3-2) peak lies in the 
nucleus, but this is not true for CO(l-O). This is very similar to another Seyfert 2 galaxy 
M51, suggesting that the AGN is playing an important role in the different behaviors of 
these two CO transition lines. 

In the spiral arms, the CO(3-2)/CO(l-0) integrated intensity ratio is well correlated 
with the star formation rate surface density, indicating that the physical conditions of 
molecular gas are related to star formation. Both the CO(3-2)/(l-0) ratio and the star 
formation rate decrease with radius from the nucleus. 
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